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Abstract 

The impact of radiative processes on kinetic equihbration is studied via a ra- 
diative transport model. The 2 o 3 processes can significantly increase the 
level of thermalization. These processes lead to an approximate coupling con- 
stant scaling of the evolution of the pressure anisotropy quahtatively differ- 
ent from the case with only 2 — > 2 partonic processes. Furthermore, thermal 
and Color Glass Condensate motivated initial conditions are shown to share 
the same asymptotic evolution when 2 -H- 3 processes are included. This 
emphasizes the unique role of radiative processes in Quark-Gluon Plasma 
thermalization. 
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1. Introduction 

The Quark-Gluon Plasma can be produced in central relativistic heavy 
ion collisions It leads to strong collective fiow and jet quenching. Ideal 

hydrodynamics is very successful in describing experimental observables in 
the low transverse momentum region |i5h8|]. Transport models can also give 
good descriptions of global observables in these collisions j9-16|. Ideal hy- 



drodynamics assumes local thermal equilibrium and the ideal hydrodynamic 



equations are valid when there is local isotropization [17(|. It is interest 



ing to study the effects of viscosity [18H24] and how the partons produced 



in heavy ion collisions thermalize |25H28|. The equilibration process can 



be studied with a microscopic transport model |29|-|32|. In particular, Xu 
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and Greiner studied the transport rates [29| and showed that radiative pro- 
cesses are important for parton momentum isotropization. A related topic is 
pressure isotropization. The Frankfurt group extended the viscous hydrody- 
namic equation to the third order |32|] for the pressure isotropization study. 
They compared pressure anisotropy evolutions starting from an isotropic ini- 
tial condition with both the extended viscous hydrodynamics and an elastic 
parton cascade with time dependent cross sections. The comparison demon- 
strates the importance of higher order corrections in describing highly viscous 
matters. In the following, we will study whether and how pressure isotropiza- 
tion depends on the inclusion of radiative processes. The pressure anisotropy 
evolutions with different initial anisotropics, energy densities, fugacities, mo- 
mentum distributions, coupling constants will be compared and contrasted 
to demonstrate the effects of radiative processes. We will further discuss how 
the system loses memory of the initial pressure anisotropy to approach the 
common asymptotic evolution and the interplay between chemical equilibra- 
tion and kinetic equilibration. 



2. Pressure Isotropization and radiative transport 

Relativistic heavy ion collisions produce hot and dense matter. The high- 
est energy density is achieved in the central cell in central collisions. Ki- 
netic equilibration in the central cell can be characterized by the pressure 
anisotropy parameter, i.e., the longitude pressure to the transverse pressure 
ratio Pl/Pt- When there is thermal equilibrium, the pressure anisotropy 
equals one. Any pressure anisotropy value that is not equal to one indicates 
non-equilibrium conditions. 

We will study the early stage when longitudinal expansion dominates 
the central cell evolution. The initial conditions will be taken to be similar 
to those expected in central relativistic heavy ion collisions. Only massless 
gluons are included and the formation proper time is set to be 0.5 fm/c. 
These gluons will be distributed uniformly inside a transverse circle with a 
radius of 5 fm and a space-time rapidity region from -5 to +5. The initial 
space-time rapidity density dN/drjg = 1000. The partonic system can start 
with an initial energy density eo = 38.2 GeV/fm^ or 76.4 GeV/fm^. In the 
local rest frame, the initial particle distribution can be isotropic or transverse. 
Both have exponential momentum distributions. In the isotropic case, eo = 
38.2 GeV/fm^ is equivalent to an initial temperature of Tq = 0.5 GeV and 
eo = 76.4 GeV/fm^ corresponds to an initial temperature of 1 GeV. The 
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transverse initial conditions are isotropic in the transverse plane with an 
initial temperature 50% higher than the isotropic case to have the same 
initial energy density. 

To study the effects of radiative processes, we will compare results with 
only the 2 — )■ 2 process that has 2 incoming gluons and 2 outgoing gluons 
with those including additional 2^3 processes. The 2—7-2 cross section is 
set to be the perturbative QCD cross section regulated by a Debye screening 
mass, i.e., 

- - ^- w 

In the above formula, as is the strong interaction coupling constant and 
fi is the Debye screening mass. The Debye screening mass is calculated 
dynamically via 

i ^i 

where V is the volume of the cell and the sum goes over all gluons in the cell. 

The effect of dynamical screening can be demonstrated by comparing 
the collision rate with the expansion rate. In a fixed box, evolution from 
a pressure anisotropy different from 1 is characterized by the collision rate 
only. However, when the system undergoes longitudinal expansion, there is 
a tendency of evolving toward pressure anisotropy. The initial pressure 
anisotropy and the expansion rate are responsible for how fast this hap- 
pens. It is the competition between collision and expansion that determines 
whether pressure isotropization can happen. For the simple case when the 
system can be approximately described by a temperature T, the cross section 
a (X l/^j? oc l/T^. The collision rate Rc = na ocT^x ~ T oc In 
the above formula, n stands for the particle number density and the relation 
between T and the proper time r for the Bjorken expansion is used. The 
volume expansion rate oc 1/t. It decreases faster than the collision rate. 
Therefore, even if the initial expansion rate is large and there is a large initial 
pressure anisotropy, resulting in a decrease of pressure anisotropy toward 0, 
the collision rate will eventually take over and pressure anisotropy will evolve 
toward 1. 

The 2^3 radiative processes are implemented for the study of the 
importance of radiative processes. In particular, the 2—7-3 cross section 
is taken to be 1/2 of the 2—7-2 cross section. This is in line with a more 
elaborate study by Xu and Greiner 15|]. To ensure the correct chemical 
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equilibrium, the detailed balance relation needs to be enforced for the 3 — )■ 2 
process. We will take all outgoing particles to be isotropic in the center of 
mass frame of the collision. This is expected to reflect the gross features of 
particle production in the central region. The 3 — )■ 2 collision rate is related 
to the reaction integral I^2i i-e., the integral over the phase space volume 
(with proper summation and averaging over internal degrees of freedom) 
of the modulus squared of the transition matrix element M. In this case, 
h2 = |i (2^)3x2 |Mp(27r)^-R2(g^'^^, 0, 0), where d = 16 is the gluon degeneracy 



factor, R2 is the two-body phase space integral as defined in [33j, and s is 
the center of mass energy squared. It is directly proportional to the 2 — )■ 3 
cross section, i.e., 

1927r2 2 
h2 = — 2 — ^23 = 127r 0-23. (3) 

The inclusion of radiative processes couples the chemical equilibration to 
the kinetic equilibration. If the system is not far from chemical equilibrium, 
then the inclusion of radiative processes will lead to additional collisions and 
isotropization. If the system is far below chemical equilibrium, there will be 
significant particle production. This will lead to smaller cross sections which 
will limit the additional isotropization relative to the case with only elastic 
processes. If the system is far above chemical equilibrium, the decrease in 
particle number will lead to larger cross sections and enhanced isotropization. 

Pressure anisotropy evolutions with elastic collisions only and with radia- 
tive processes are compared in Fig. [H First, let us focus on panel (b) and 
discuss some of the general features of the evolutions. There are two curves 
for the case including 2^3. One is for the case starting with an initial 
isotropic distribution, i.e., initial Pl/Pt = 1- The other is for the case with 
an inside-outside type of initial distribution, i.e., initial Pl/Pt = 0. Even 
though they start from quite different initial pressure anisotropies, they ap- 
proach the same asymptotic evolution at late times. In other words, the 
memory of initial anisotropy is lost after some characteristic time. The case 
with only 2 — )■ 2 collisions has approximately the same behavior. The two 
curves approach and meet each other at late times. This feature is the same 
for all other cases. The importance of radiative processes is reflected in the 
difference of the radiative case and the case with only elastic collisions. The 
radiative case has a pressure anisotropy of about 0.78 at 2 fm/c, significantly 
larger than that of the elastic only case of 0.65. This shows that radiative 
processes can significantly enhance thermalization. 
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Figure 1: Pressure anisotropy evolution. Solid curves are for the case including the 2 o 3 
processes while dashed curves are for the case with 2^-2 only. The upper panels have 
initial energy density eo — 38.2 GeV/fm"^ while the lower panels have eo = 76.4 GeV/fm"^. 
The left panels are calculated with Os = 0.3 while the right panels have as — 0.6. 
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When as decreases from 0.6 to 0.3, the asymptotic value at 2 fm/c de- 
creases for both the case with 2 3 and without 2 4-7-3. When the initial 
energy density increases from 38.2 GeV/fm^, the case with 2 -f-)- 3 increases 
slightly leading to an approximate as scaling insensitive to the initial en- 
ergy density. In contrast, the case with only the 2 — )■ 2 process increases 
drastically. Because of the slight increase in the radiative case and drastic 
increase in the elastic only case as the energy density increases, when the 
initial energy density eo = 76.4 GeV/fm^, the radiative and elastic only cases 
have about the same pressure anisotropy at 2 fm/c. The elastic only case is 
sensitive to both the initial energy density and the coupling constant. The 
large initial energy density and small coupling constant case is related to 
the small initial energy density and large coupling constant case via the eoa^ 
scaling. The elastic only curves in panel (b) and panel (c) follow this eoa^ 
scaling and are identical for these two cases. 

The above analysis demonstrates that evolutions with and without radia- 
tive processes have different dependences on system parameters. In the case 
with radiative processes, when the energy density is increased, the screening 
mass decreases, increasing the cross section and collision rate. However, this 
is counteracted by the additional particle production. As the system falls far 
below chemical equilibrium, additional particle production tends to increase 
the screening mass and reduce the collision rate. Lacking this balancing fac- 
tor, the elastic only case cannot maintain approximately the same evolution 
as with the lower energy density case and the pressure anisotropy increases 
drastically. Decreasing the coupling constant can compensate for the increase 
in pressure anisotropy caused by the energy density increase. As the initial 
binary collision cross section is proportional to e^as with fixed initial num- 
ber density, when the coupling constant is decreased by the same factor as 
the energy density increase, there is an exact e^ag scaling in the elastic only 
case. This is qualitatively different from the approximate as scaling in the 
radiative case. 

The two data sets with different initial pressure anisotropics but otherwise 
same parameters seem to converge toward a common final evolution. The 
difference between these two data sets can be used to study the rate of 
convergence as shown in Fig. |2]for cases including the 2 <H- 3 processes. The 
difference quickly decreases and at late times fits well with an exponential 
decrease with proper time. In the case with 2 — )■ 2 only, the two curves may 
cross each other at late times (see e.g. panel (d) in Fig. [1]). In any case, by 
2 fm/c, the differences all go to about or below 1% of the initial difference. 
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T (fm/c) 

Figure 2: Pressure anisotropy difference evolution. These results all include the radiative 
processes. Crosses (squares) are for initial energy density cq — 38.2(76.4) GeV/fm^. Solid 
and dashed curves are exponential fits for r > 0.8 fm/c. For the two data sets with the 
same symbol, the upper one is for — 0.3 and the lower one is for as = 0.6. 
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Figure 3; Proper time evolution of fugacity A for initial space-time rapidity density 
dNo/drjs — 1000 and coupling constant as = 0.6. Lines are for isotropic exponential 
initial conditions and points are for transverse exponential initial conditions. Solid lines 
and circles are results including 2 o 3. Dashed lines and pluses are results with 2^-2 
only. For the two data sets with the same symbol, the upper one is for cq = 38.2 GeV/fm'^ 
and the lower one is for eq = TGA GeV/fm^. 

The system approaches the asymptotic evolution on a very short time scale. 
The comparison in Fig. |2] shows that the approach to asymptotic evolution 
is faster with higher initial energy density or higher coupling constant. It is 
interesting to notice that within error bars, these lines seem to go through 
the same point. We also looked at the curves with only elastic processes. 
They also seem to go through a common point, even though that point is 
different from the one with inelastic processes. 

To get a better understanding of the relation between the observed pres- 
sure isotropization and radiative processes, it is helpful to look at the fugacity 
evolution. Define fugacity to be the ratio of particle density to the equilib- 
rium particle density, i.e., A = n/neq- When the temperature T is taken 
to be the average kinetic energy per degree of freedom, i.e., T = e/(3n), 
the resultant Ugq leads to the expression A = 277r^?T,^/(16e^). For systems 
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in thermodynamical equilibrium, i.e., in chemical and thermal equilibrium, 
A = 1. The comparison of fugacity evolutions in Fig. E] demonstrates the in- 
terplay between thermal equilibration and chemical equilibration. Transverse 
initial conditions have evolutions lower than their isotropic initial condition 
counterparts mainly because of their slower early energy density evolutions 
caused by the inside-outside space-momentum correlation. When only 2—7-2 
is included, the system falls farther and farther away from chemical equilib- 
rium. As the 2 -H- 3 processes are turned on, the system is able to approach 
chemical equilibrium. Even for transverse initial conditions with eo = 38.2 
GeV/fm^ where there is a decrease in fugacity mainly due to longitudinal 
expansion, at late times, the fugacity is able to reach a value that is com- 
parable to other cases. Note that in this case, there are more annihilations 
than productions initially even when A ~ 1 because of the transverse spa- 
tial distribution. This is why the early evolution is slightly lower than the 
corresponding elastic only case. It demonstrates that chemical equilibration 
depends on kinetic equilibration in the radiative case. For various cases with 
the 2 o 3 processes, the final fugacities at 2 fm/c are all on the 80% to 90% 
level. Note that the eo = 38.2 GeV/fm^ case has an initial fugacity close to 
1, while the eo = 76.4 GeV/fm^ case is significantly undersaturated initially. 
As discussed before, this difference in initial fugacities is behind the observed 
different energy dependences of pressure anisotropy evolutions for the elastic 
only and the radiative cases. 

3. Summary and discussions 

The above study demonstrates the different behaviors of evolutions with 
and without radiative processes. When radiative processes are included, 
there is an approximate scaling that is in contrast with the eoa^ scal- 
ing seen when only elastic processes are included. Evolutions with different 
initial pressure anisotropics appear to approach the same asymptotic evo- 
lution exponentially at late times. If a partonic system is initially close to 
chemical equilibrium, radiative processes can significantly enhance pressure 
isotropization. For systems significantly undersaturated at the initial time, 
evolutions with and without radiative processes can be close to each other. 
This dependence on initial fugacity indicates that it may be necessary to use 
different K factors for elastic simulations at different collision energies. 

The above calculations start with initial exponential momentum distribu- 
tions. The dependence on the initial momentum distribution can be studied 



9 



1.2 
1 

0.8 

9:" 0.6 
Q- 

0.4 
0.2 


0.5 1 1.5 2 
T (fm/c) 

Figure 4: Pressure anisotropy evolution for dNa/drfs — 1000, eq = 38.2 GeV/fm'^, = 0.6. 
Lines are for exponential initial momentum distributions and points are for condensate 
initial momentum distributions. Solid lines and circles are for the case including the 2 -o- 3 
processes while dashed lines and pluses are for the case with 2 — s> 2 only. 
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by comparing the results with those that start with initial momentum dis- 
tributions motivated by the Color Glass Condensate. Fig. H] has additional 
calculations with step function initial momentum distributions. A clear en- 
hancement of pressure isotropization is observed with the step function initial 
distributions when only elastic processes are included. When inelastic pro- 
cesses are included, there is not much change in the pressure anisotropy evo- 
lution. In other words, the pressure isotropization is robust against changes 
in the initial momentum distribution. This difference between the elastic 
only and with radiative reflects the faster thermalization when radiative pro- 
cesses are included. Without them, the step function momentum distribution 
maintains its shape for a longer period of time, resulting in a smaller screen- 
ing mass and larger collision rate per particle and a significant change in 
the pressure anisotropy evolution. This comparison demonstrates again the 
importance of radiative processes in microscopic simulations of relativistic 
heavy ion collisions. 
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